Kinetics of photoinduced ordering in azo-dye films: 
two-state and diffusion models 



Alexei D. Kiselev,^'^'B Vladimir G. Chigrinov/'[l| and Hoi-Sing Kwok^'H 

^ Hong Kong University of Science and Technology, 
Clear Water Bay, Kowloon, Hong Kong 
^Institute of Physics of National Academy of Sciences of Ukraine, 
prospekt Nauki 46, 03028 Kyi'v, Ukraine 
(Dated: December 10, 2008) 

We theoretically study the kinetics of photoinduced ordering in azo-dye photoalign- 
ing layers and present the results of modeling performed using two different phe- 
nomenological approaches. A phenomenological two state model is deduced from 
the master equation for the one-particle distribution functions of an ensemble of 
two-level molecular systems by specifying the angular redistribution probabilities 
and by expressing the order parameter correlation functions in terms of the order 
parameter tensor. Using an alternative approach that describes light induced reori- 
entation of azo-dye molecules in terms of a rotational Brownian motion, we formulate 
the two-dimensional (2D) diffusion model as the free energy Fokker-Planck equation 
simplified for the limiting regime of purely in-plane reorientation. The models are 
employed to interpret the irradiation time dependence of the absorption order pa- 
rameters defined in terms of the the principal extinction (absorption) coefficients. 
Using the exact solution to the light transmission problem for a biaxially anisotropic 
absorbing layer, these coefficients are extracted from the absorbance-vs-incidence 
angle curves measured at different irradiation doses for the probe light linearly po- 
larized parallel and perpendicular to the plane of incidence. It is found that, in 
the azo-dye films, the transient photoinduced structures are biaxially anisotropic 
whereas the photosteady and the initial states are uniaxial. 

PACS numbers: 61.30.Gd, 42.70.Gi, 82.50.Hp 
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I. INTRODUCTION 

It has long been known that some photosensitive materials such as compounds containing 
azobenzene and its derivatives may become dichroic and birefringent under the action of 
light. This phenomenon — the so-called effect of photoinduced optical anisotropy (POA) — 
has a long history dating back almost nine decades to the paper by Weigert [1!]. 

The Weigert effect (POA) has been attracted much attention over the past few decades 
because of its technological importance in providing tools to produce the light-controlled 
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anisotropy. For example, the materials that exhibit POA are very promising for use in 
many photonic applications jl, 0, 0, 0] • 

It is also well known that producing substrates with anisotropic anchoring properties is 
one of the key procedures in the fabrication of liquid crystal electrooptic devices. The tra- 
ditional method widely used to align liquid crystal display cells involves mechanical rubbing 
of aligning layers and has a number of the well known difficulties 0. The photoalignment 
technique suggested in Refs. is an alternative method that avoids the drawbacks of 

the mechanical surface treatment by using linearly polarized ultraviolet (UV) light to induce 



anisotropy of the angular distribution of molecules in a photosensitive film [10|, . Thus 
the phenomenon of POA (the Weigert effect) is at the heart of the photoalignment method. 

Light induced ordering in photosensitive materials, though not being understood very 
well, can generally occur by a variety of photochemically induced processes. These typically 
may involve such transformations as photoisomerization, crosslinking, photodimerization 



and photodecomposition (a recent review can be found in Ref. [Uj, Il2l]). 

So, the mechanism underlying POA and its properties cannot be universal. Rather they 
crucially depend on the material in question and on a number of additional factors such as 
irradiation conditions, surface interactions etc. In particular, these factors combined with 
the action of light may result in different regimes of the photoinduced ordering kinetics 
leading to the formation of various photoinduced orientational structures (uniaxial, biaxial, 
splayed) . 

POA was initially studied in viscous solutions of azodyes [l3| and in azodye-polymer 
blends where the anisotropy was found to be rather unstable. This is the case where the 
photoinduced anisotropy disappears after switching off the irradiation 13, 14, H, 16, 13, 18 . 
By contrast to this case, POA can be long term stable. 

The stable POA was observed in polymers containing chemically linked azochromophores 
(azopolymers) ^]. It turned out that stable anis otropy can be induced in both amorphous 



and liquid crystalline (LC) azopolymers jl, Is], 111, 20. l2ll. [2I . 23 

The photoalignment has also been studied in a number of similar polymer systems in- 
cluding dye doped polymer layers 0:_24|, cinnamate polymer derivatives [s], 0, 



26 and 



side chain azopolymers I, |22|, |27j. In addition, the films containing photochemically 



stable azo dye structures (azobenzene sulfuric dyes) were recently investigated as new pho- 



29, 30 



toaligning materials for nematic liquid crystal (NLC) cells [2 

In Ref. [11], it was found that, owing to high degree of the photoinduced ordering, these 
films used as aligning substrates are characterized by the anchoring energy strengths com- 
parable to the rubbed polyimide films. For these materials, the voltage holding ratio and 
thermal stability of the alignment turned out to be high. The azo-dye films are thus promis- 
ing materials for applications in liquid crystal devices. 

According to Ref. [13], the anchoring characteristics of the azo-dye films such as the polar 
and azimuthal anchoring energies are strongly infiuenced by the photoinduced ordering. In 
this paper the kinetics of such ordering will be of our primary interest. More specifically, 
we deal with theoretical approaches and related phenomenological models describing how 
amount of the photoinduced anisotropy characterized by absorption dichroism evolves in 
time upon illumination and after switching it off. 
There are a number of models 



17, 23 



33, M 35 



exhibiting POA driven by the trans- cis photoisomerization. 
sample is treated as an ensemble of two level molecular systems 
characterized by elongated rod-like molecular conformation can be regarded as the ground 



formulated for azocompounds 
Generally, in these models, a 
the stable trans isomers 



state molecules whereas the bent banana-like shaped cis isomers are represented by the 
excited molecules. 

The photoisomerization mechanism assumes that the key processes behind the orienta- 
tional ordering of azo-dye molecules are photochemically induced trans- cis isomerization 
and subsequent thermal and/or photochemical cis-trans back isomerization of azobenzene 
chromophores. 

Owing to pronounced absorption dichroism of photoactive groups, the rate of the photoin- 
duced isomerization strongly depends on orientation of the azo-dye molecules relative to the 
polarization vector of the actinic light, E[/y. Since the optical transition dipole moment is 
approximately directed along the long molecular axis, the molecules oriented perpendicular 
to E[7y are almost inactive. 

When the cis isomers are short-living, the cis state becomes temporary populated during 
photoisomerization but reacts immediately back to the stable trans isomeric form. The 
trans- cis-trans isomerization cycles are accompanied by rotations of the azo-dye molecules 
that tend to minimize the absorption and become oriented along directions normal to the 
polarization vector of the exciting light Euv- Non- photoactive groups may then undergo 
reorientation due to cooperative motion |17|, [l9|, [33^ 



The above scenario, initially suggested in Ref. [l3j, is known as the regime of photoorien 



tation (angular redistribution) where the lifetime of cis isomers is short and POA is mainly 
due to the angular redistribution of the long axes of the trans molecules during the trans-cis- 
trans photoisomerization cycles. Note that, in the opposite case of long-living cis isomers, 
the regime of angular hole burning (photoselection) occurs so that the anisotropy is caused 
by angular selective burnin g o f mesogenic trans isomers due to stimulated transitions to 



non-mesogenic cis form [22l, l35|, l37 



From the above it might be concluded that, whichever regime of the ordering takes 
place, the photoinduced orientational structure results from preferential alignment of azo- 
dye molecules along the directions perpendicular to the polarization vector of the actinic 
light, E[/y, determined by the dependence of the photoisomerization rate on the angle 
between E[/y and the long molecular axis. So, it can be expected that the structure will be 
uniaxially anisotropic with the optical axis directed along the polarization vector. 

Experimentally, this is, however, not the case. For example, constraints imposed by 
a medium may suppress out-of-plane reorientation of the azobezene chromophores giving 
rise to the structures with strongly preferred in-plane alignment [l^]. Another symmetry 
breaking effect induced by polymeric environment is that the photoinduced orientational 



structures can be biaxial [2l|, |23|, |35|, l38|, |39|, |40(| (a recent review concerning medium effects 



on photochemical processes can be found in |4l[|). 

It was recently found that, similar to the polymer systems, the long-term stable POA in 
the azo-dye SDl films is characterized by the biaxial photoinduced structures with favored 
in-plane alignment Unlike azopolymers, photochromism in these films is extremely 

weak so that it is very difficult to unambiguously detect the presence of a noticeable fraction 
of cis isomers. 

As compared to the polymer systems, modeling of photoinduced ordering in the azo-dye 
films has received little attention. In this paper we intend to fill the gap and describe the 
symmetry breaking and biaxiality effects using phenomenological models formulated on the 
basis of a unified approach to the kinetics of POA j23|, 3^]. The layout of the paper is as 
follows. 

In Sec. Ill Al we introduce necessary notations and discuss the relationship between the 



order parameter and the absorption tensors. Then, in Sec. Ill B\ we recapitulate the the- 
ory '35| by assuming that the azo-dye molecules can be represented by two level molecular 
systems. This theoretical approach is based on the master equation combined with the ki- 
netic equation for the additional (matrix) system, which phenomenologically accounts for 
the presence of long-living anisotropic (angular) correlations. 

In Sec. Ill CI a phenomenological two state model is introduced by specifying the angular 
redistribution probabilities and by expressing the order parameter correlation functions in 
terms of the order parameter tensor. In this model, the regime of photoorientation with 
short living excited molecules is characterized by weak photochromism and negligibly small 
fraction of cis isomers that rapidly decays after switching off irradiation. 



According to Ref. [43|, when the photochemical processes underlying photoisomerization 
are hindered, the process of photoinduced reorientation can be alternatively described as 
rotational diffusion of azo-dye molecules under the action of the polarized light. 

In Sec. IIII At we show that diffusion models of POA can be formulated as the free energy 
Fokker-Planck equation [43] describing light induced reorientation of azo-dye molecules as 
rotational Brownian motion governed by the effective mean field potential. Using this ap- 
proach, the diffusion model suggested in j43| can be easily extended to the case of biaxial 
orientational structures. In Sec. IIII B\ we introduce and study the simphfied two-dimensional 
(2D) diffusion model that can be regarded as the first approximation representing the regime 
of purely in-plane reorientation. 

The two state and 2D diffusion models are employed to interpret the experimental data 
in Sec. IIVI Finally, in Sec. |V] we present our results and make some concluding remarks. 
Technical details on solving the light transmission problem for a biaxially anisotropic ab- 
sorbing layer and on using the analytical result to extract the extinction coefficients from 
the measured dependence of absorbance on the incidence angle are relegated to Appendix. 



II. MASTER EQUATION AND TWO-STATE MODELS 

A. Order parameters, absorption tensor and biaxiality 

We assume that azo-dye molecules are cylindrically symmetric and orientation 
of a molecule in the azo-dye film can be specified by the unit vector, u = 
(sin ^ cos 0, sin 6* sin 0, cos 6*), directed along the long molecular axis. Quadrupolar orienta- 
tional ordering of the molecules is then characterized using the traceless symmetric second- 



rank tensor 45 



Q(u) = (3u®u-I)/2, (1) 

where I is the identity matrix. 

The dyadic ([T]) averaged over orientation of molecules with the one-particle distribution 
function p(r,u), describing the orient at ion- density profile of azo-dye molecules, is propor- 
tional to the order parameter tensor S(r) 



/ 



p(r,u)Q(u)du = p(r)S(r), (2) 



where du = sin^d^d^, p(r, u) = p(r)/(r, u), p(r) = J p(r, u)du is the density profile and 
/(r, u) is the normalized angular distribution. 




FIG. 1: Frame of reference: the z axis is normal to the substrate and the polarization 
vector of the activating light is directed along the y axis. 



Throughout the paper we restrict ourselves to the case of spatially homogeneous systems 
with p(r, u) = p/(u). For such systems, the order parameter tensor is given by 



(Q) = 1 Q(u)/(u) du = S, 

S = 5(3 d eg) d - I)/2 + P(m ® rh - 1 ® i)/2, 



(3) 
(4) 



where Ai = 5", A2 = — (-P + >S')/2 and A3 = — (Ai + A2) = (P — 5')/2 are the eigenvalues of 
the order parameter tensor S; the eigenvector d corresponding to the largest in magnitude 
eigenvalue, Ai = S", \S\ = max{|Ai|, IA2I, IA3I}, is the unit vector known as the director; P 
is the biaxiality parameter and the eigenvectors {d, m, 1} form a right-handed orthonormal 
tripod. 

In our case, the y axis is directed along the polarization vector of the activating UV light, 
E = Ey, the z axis is normal to the substrates and the unit vector x = [y x z] is parallel 
to the X axis (see Fig. [T]). On symmetry grounds, it can be expected that the basis vectors 
{x, y, z} define the principal axes of the order parameter tensor. So, the tensor is given by 



S = diag(5'^, Sy, 5^) = 5^. x (g) x 



(5) 



Then the dielectric tensor, e, can also be written in the diagonal form 



£ — dia.g{ex,ey,ez), Sa/3 — ^a^ap- 



(6) 



In the presence of absorption, the tensor is complex-valued and its principal values, 
{e^,., ey, e^}, are expressed in terms of the refractive indices, {nx \ ny \ riz^} and the extinction 
coefficients, {k^,^^,^^}, as follows [46| : 



+ ie'^i fiea = = {n'^a^ + ii^af- (7) 



We can now define the absorption order parameters through the relation 



S) = —, T = -r\ 7-^ 7^ ? 7^ (8) 

2{k, + + Kk) 2[Dt^ + + Df -' ' ^ ^ ^ ^ ^ 



where the optical densities {Z^i"'*, D^'^'', Di"^} are proportional to the extinction coefficients: 
D-"'' oc Ki. Note that the optical density D^^^ = Dy""^ [D^^^ = Z^i""*] can be determined 
experimentally by measuring the absorption coefficient for a testing beam which is propa- 
gating along the normal to the film substrate (the z axis) and is linearly polarized parallel 
[perpendicular] to the polarization vector of the activating UV light (the y axis). 



Now, following Ref. [30|, |40|, we dwell briefly on the relation between the orientational 
and the absorption order parameters deflned in Eq. ([5]) and Eq. ([H]), respectively. To this 
end, we begin with the absorption tensor of an azo-dye molecule 

o-ij{u) = cr±6ij + (a-|| - a±)uiUj, (9) 

which is assumed to be uniaxially anisotropic. Its orientational average takes the following 
matrix form 

(cr) = (a,,I + 2AaS)/3, (10) 
o"av = + 2(7_L, Aa = a\\-a±, (11) 

where the angular brackets (. . .) denote orientational averaging (see Eq. ([3])). 

In the low concentration approximation, the optical densities are proportional to the 
corresponding components of the tensor (fTOjl 



d\^^ =D^:^ oc p(aav + 2 Aa5,.)/3, (12a) 

=4'^) oc p((Tav + 2 Aa5,)/3, (12b) 
ocp(aav + 2A(T5,)/3 (12c) 

and on substituting the expressions for the optical densities f lT^ into Eq. ([H]) we obtain 

St'^=r^S,, (13) 

where = Aa/dav = cTq/ (3 + da). So, the absorption order parameters ^ are equal to the 
corresponding elements of the order parameter tensor ^ only in the limiting case where 
absorption of waves propagating along the long molecular axis is negligibly small: — 
and (Jav = 3cr_L + Aa — > Act. Note that the average optical density D^x'' + D^"^ + Z^i'*^ is 
proportional to paav and thus typically does not depend on the irradiation dose. 



B. Master equation 



We shall assume that the azo-dye molecules can be represented by the two-level molecular 
systems with the two states: the ground state and the excited state. Angular distribution 
of the molecules in the ground state at time t is characterized by the number distribution 
function NG{u,t) = VpG{u,t), where V is the volume and Pg(u, t) is the corresponding 
one-particle distribution function. 

Similarly, the azo-dye molecules in the excited state are characterized by the function: 
iVE(u, t) = VpE{u,t). Then the number of molecules in the ground and excited states is 
given by 



Nait) = Nncit) = j Ng{vl, t) du = V j pg(u, t) du, 



(14) 



NE{t) = NuEit) = j Ne{u, t)du = V j pe(u, t) du, ndt) + nE{t) = 1, (15) 

where is the total number of molecules; and Ue are the concentrations of non-excited 

du = / d(f) sin 9 d9 . 
Jo Jo 

The normalized angular distribution functions, /^(n, t), of the ground state {a = G) and 
the excited {a = E) molecules can be conveniently defined through the relation 



(16) 



linking the one-particle distribution function, pa, and the corresponding concentration, n^. 

The presence of long-living angular correlations coming from anisotropic interactions 
between azo-dye molecules and collective modes of confining environment can be taken 

I47I . In this 



into account by using the phenomenological approach suggested in Refs. |4( 
approach, the effective anisotropic field, that results in the long-term stability effect and 
determines angular distribution of the molecules in the stationary regime, is introduced 
through the additional angular distribution function, /^(u, t). characterizing the additional 
subsystem that, for brevity, will be referred to as the matrix system. 

It bears close resemblance to the equilibrium distribution of the mean field theories of 
photoinduced optical anisotropy [3l|, 32, HI]. In these theories, this distribution has been 
assumed to be proportional to exp(— \^(u)//cbT), where V{u) is the mean-field potential 
that depends on the order parameter tensor. 

We shall write the kinetic rate equations for Na{u,t) in the general form of master 



equation [48|, |49|, |50 



dt 

+ 7o 



dN^ 
dt 



+ 5Z / [ ^ I ^i')^/3(u', t) - W{p, u I a, u) N^{u, t) 



Diff 

Nna{t) / r„_„(u, u )/m(u ,t)du' - iV„(u,)f:) 



du' 



(17) 



where a,/? G {G, E}. 

The first term on the right hand side of Eq. ( ITTl) is due to rotational diffusion of azo-dye 
molecules. In frictionless models this term is absent. It will be considered later on in Sec. Illli 



Now we need to specify the rate of the G — > E transition stimulated by the incident UV 
hght. For the electromagnetic wave linearly polarized along the y-axis the transition rate 



can be written as follows [16|, l51 



l^(E, u I G, u') = re-,(u, u') Pg(u'), (18) 

P,(U) = {nWt)-'<!>g^eJ2''if(^)^^^J=l9luv{l + (TaUl) 

id 

= qgIuv{3 + (Ta + 2aaQyy{u))/3 (19) 

where a^^\u) is the tensor of absorption cross section for the molecule in the ground state 
oriented along u: a^f = cr^f'Sij + (cr||^'' — a^f^) UiUj; a a = (cX||^^ — a'f^)/a^f' is the absorption 
anisotropy parameter; hjjt is the photon energy; ^g^e is the quantum yield of the process 
and rc,_e(u, u') describes the angular redistribution of the molecules in the excited state; / 

is the pumping intensity and qg = {f}Wt)~^^g^eO'^±' ■ 

Similar line of reasoning applies to the E — G transition to yield the expression for the 
rate: 

W{G, u I E, u') = TeF^i (u, u') + qJuv T^"? (u, u') , (20) 

where qe = {}TWt)^^^e~>gO'^'^'' and 7e = I/ts, is the lifetime of the excited state and the 
anisotropic part of the absorption cross section is disregarded, cri^i*^^ = a^f^ = a^^\ 

Equation (l20l) ) implies that the process of angular redistribution for spontaneous and 
stimulated transitions can be different. All the angular redistribution probabilities are nor- 
malized so as to meet the standard normalization condition for probability densities: 

r^_,(u,u')du= 1. (21) 

Using the system ffT7|) and the relations ffT^ - fl^Uj) it is not difficult to deduce the equation 
for nG{t)\ 

= 7e (1 - ^g) - {Pg)GnG , 7e = 7e + qehv , (22) 

where the angular brackets (...)« stand for averaging over the angles with the distribution 
function /„ . Owing to the condition ( 12T1) . this equation does not depend on the form of the 
angular redistribution probabilities. 

The last square bracketed term on the right hand side of (ITTI) describes the process that 
equilibrates the absorbing molecules and the matrix system in the absence of irradiation. 
The angular redistribution probabilities rQ,_m(u, u') meet the normalization condition, so 
that thermal relaxation does not change the total fractions A^^g and A^e- If there is no angular 
redistribution, then rQ,_m(u, u') = 5(u — u') and both equilibrium angular distributions /^''^ 
and f^''^ are equal to /„• 

The latter is the case for the mean field models considered in 31, 3^, 34 . In these mod- 



els the excited molecules [cis fragments) are assumed to be long-living with 7e = and 
7g = 7e- We can now recover the models by setting the angular redistribution probabilities 
rg_e(u, u') and re_g(u, u') equal to the equilibrium distribution, = p(u), determined by 
the mean-field potential ^(u): p(u) oc exp(— l^/fc^T). So, the mean field approach intro- 
duces the angular redistribution operators acting as projectors onto the angular distribution 



of the matrix system. This is the order parameter dependent distribution that characterizes 
orientation of the azo-molecules after excitation. 

An alternative and a more general approach is to determine the distribution function 
/m(u, t) from the kinetic equation that can be written in the following form fSSj: 



dfm{ii,t) 
dt 



Q={G,E} 



(23) 



Equations for the angular distribution functions /g(u, t) and /e(u, t) can be derived 
from ffTTl) by using the relations (fTS|l - fl^ . The result is as follows 



dt 
- 7e 



-HQ 



Pg)Gh - / r,_,(u,u')Pg(u')/G(fi',t)du' 



/e - / re-m(u, u')/m(u', t) du' 



(24) 



dt 



-TiG 



[P,(u)-(P,)G]/G+7eriE j r(^_^l(u,u')/E(u',t)du' 



- (7e + Qel) riE fc + Qeln^ I r^?f^(u, u')/e(u', t) du' 



7g ^g 



/g - / rg_^(u,u')/m(u',t) du' 



(25) 



The system of equations (1^ and can be used as a starting point to formulate 

a number of phenomenological models of POA. We have already shown how the mean 
field theories of jU, 32, can be reformulated in terms of the angular redistribution 
probabilities. 



C. Two-state model 

We can now describe our two state model. To this end, we follow the line of reasoning 



presented in Refs. [351, |40 



In this model, the angular redistribution probabilities T^-g and F^-e are both assumed 
to be isotropic: 

r,_,(u, u') = r(^_^)(u, u') = vfJX^, = i = f-- (26) 

Since we have neglected anisotropy of the excited molecules, it is reasonable to suppose that 
the equilibrium distribution of such molecules is also isotropic, /^'^'^^ = fiso = (4vr)~^, so that 

7E = 7^"^^ = 0. (27) 
From the other hand, we assume that there is no angular redistribution 

rc,_m(u, u') = r„_c,(u, u') = 5(u - u'), (28) 



and the equilibrium angular distribution of molecules in the ground state is determined by 
the matrix system: f^"^^ = /„. 

Equilibrium properties of excited and ground state molecules are thus characterized by 
two different equilibrium angular distributions: ftso and fm, respectively. It means that in 
our model the anisotropic field represented by fm does not influence the angular distribution 
of non-mesogenic excited molecules. 



dS^^j-' (G) ~ \c(G), 

^G—Qf - = - — Qghv no G\j.'yy - 7e (1 - ncjSlj '+ 

+ 7GnGiSif -Sff^), + (29a) 

^ = -7^^0(5(7) -^f), (29b) 

(G) 

where G\j.\^^ is the order parameter correlation function given by 



G'SU = (g.,(u)Q™„(ii))G - 4""^ S^^l (30) 



The key point of the approach suggested in Ref. [35[ is the assumption that the cor- 
relators fl30p which characterize response of azo-dye to the pumping light and enter the 
kinetic equations for the order parameter components (I29al) . can be expressed in terms of 
the averaged order parameters Si^ . 

riG^ ~ %IuvnG (5/7 + 2A/7 S - \'S') 

-7e(l -^0)6* + 70^0(5*^ - 5), (31a) 

^G^^ = - ^ qgluvncK'^n + AS) AS 



dASm 
dt 



7e (1 - nG)AS + ^cnciASm - AS), (31b) 
ImnciSm-S), (31c) 

'jmnciASm - AS) , (31d) 



where 7m — 'yrn \ S — Sxx\ AS — Syy^ Szz \ Sm — Sxx^ and A Sm — Syy"^ Szz"^ ■ 

It was shown that the parabolic approximation used in Ref. [23| can be improved by 
rescaling the order parameter components: S^^^ ^S^^^ with A = (1 + 0.6-\/30)/7 com- 
puted from the condition that there are no fluctuations provided the molecules are perfectly 
aligned along the coordinate unit vector ej: G\fl^ = at S^^^^ = 1. In Ref. |35| this heuris- 
tic procedure has also been found to provide a reasonably accurate approximation for the 
correlators calculated by assuming that the angular distribution of molecules can be taken 
in the form of distribution functions used in the variational mean field theories of liquid 



crystals |45l. 152 



1. Long term stability and photo steady state 



Mathematically, our model is described by equations for the order parameters and the 
concentration given in Eq. ( !3T|) and Eq. ( !22l) . respectively. We may now pass on to discussing 
some of its general properties. 

Our first remark concerns the effect of the long term stability of POA. It means that there 
is the amount of the photoinduced anisotropy preserved intact for long time after switching 
off the light. Clearly, this is a memory effect and the system does not relax back to the off 
state characterized by irradiation independent equilibrium values of the order parameters. 

In order to see how this effect is described in our model, we assume that the activating 
light is switched off at time t = tofr and consider subsequent evolution of the order parameters 
at t > toff- In the absence of irradiation, Eq. (IHTI) decomposes into two decoupled identical 
systems of equations, for two pairs of the order parameters: {S, Sm} and {AS, ASm}- So, 
without the loss of generality we may restrict ourselves to the evolution of the x components 
of order parameters, {S, Sm} = {Si^\ s!cx^} governed by the equations 

(99 

— = -7e (l/nc - \)S + 7G(^n^ - S), (32) 
dS 

= ~lm{Sm - S) (33) 

supplemented with the initial conditions 

'S'(tofr) = 'S'ofT, Sm{tof[) = sl^\ (34) 

At Iifv = 0, equation for the concentration fl^ is easy to solve. So, for the initial value 
problem with nQ{tos) = rios, we have 

1 - no = exp(-7eAt)(l - nofi). At = t - tos- (35) 

From Eq. flH31) it is clear that, in the limiting case of short living excited state with 
At ^ l/7e5 the concentration of excited molecules, = I — uq, rapidly decays to zero. In 
this regime, the first term on the right hand side of Eq. ( !32l) is negligibly small and can be 
disregarded. Equations ( l36l) and ( 1371) can now be easily solved to yield the formulas 

Sit) = Sif + ^g[SoS - exp(-7At)/7, (36) 
SUt) = - 7m[5off - ^iff ^] exp(-7At)/7, (37) 

where 7 = 7m + 7g and 

sS^ = {lmSo, + loSiP)h. (38) 

Evidently, the order parameters defined in Eqs. (1361) and (1371) evolve in time approaching the 
stationary value (l38l) . The memory effect manifests itself in the dependence of the stationary 
order parameter, S^^^\ on the (initial) conditions ( IMI) at the instant the activating light is 
switched off. 

The photosteady states reached in the long irradiation time limit are represented by 
stationary solutions of the system (PUI) and the concentration equation 



The steady state concentration of excited molecules can be expressed in terms of the 
steady state order parameter, Syy \ through the relation 



G 3(r + l) + a,(l + 24f ^ ' 



where r = 'je/{qgluv)- From Eq. (I39l) it can be seen that, in the case where the life time of 
the excited state is short and the ratio r is large, r ^ 1, the fraction of the excited molecules 
is negligible, so that n^Q^ ^ 1. 

On substituting Eq. ffT^ into the steady state relation 

(Q,,(u)P,(u))g = (40) 
derived from Eqs. fl29al) and fl22l) we obtain equation for the steady state order parameters 

5ff = {Q^.{n)Qyy{n))G. (41) 

From Eq. (129bp the difference between the order parameters of the matrix system and 
the ground state azo-dye molecules dies out as the photosteady state is approached, 
S^f^ — Slj^^ — > at t — > cxD. Interestingly, Eq. (HTl) shows that the order parameters, 
Slf\ characterizing the regime of photosaturation are independent of the light intensity. 

For the specific form of the correlators used to obtain the system (!HT]) . the photosteady 
state is uniaxial with Syy^ = Sif^ and the x component of the order parameter tensor, 
'S'st = Sxx\ can be found by solving the equation 

2aa (1/5 + 2A/7 5,t - A^^^j ^ ^^^^g ^ ^^^^ ^ 25,0). (42) 



III. NONLINEAR FOKKER-PLANCK EQUATIONS AND DIFFUSION MODEL 

A. Mean field Fokker- Planck equations 



In this section we extend the diffusion model 4^ by using the approach based on Fokker- 
Planck (F-P) equations of the following general form (ssl ]: 



di 



ap = a 



-D 



(drft) p ^(diff) 



P 



a 



(43) 



d 



(drft) 



where a = ^— and P(x, t) is the probability density (distribution function); D 

drift vector and D^'j^^^ is the diffusion tensor. 

As opposed to the linear case, in nonlinear F-P equations, either the drift vector, D 
or the diffusion tensor, D^^^^^ Ar^^^^^A r^,^ ^;cfv.iK,,fir^v^ A,^^fi^v^ c n('^''^*) — n('^''^*) 



is the 



(drft) 



and Dlj™' = -D^^'^"* [Pi. The theory and applications of such equations were recently reviewed 
in the monograph [4^. Interestingly, according to Refs. 5^, l55|], nonlinear F-P equations 



, depend on the distribution function, P: _D^^""'-' = /^^^""'-'[P] 



are derived by approximating the master equation with nonhnear effects introduced through 
the generahzed transition rates. 

More specifically, we concentrate on the special case of the so-called nonlinear mean-field 
Fokker- Planck (F-P) equations 

a,p.,{p.,^^ffl}.v.{pV^fa} (44, 

that are characterized by the effective free energy functional, i^[-P]. In addition, the case of 
the rotational Brownian motion will be of our primary interest. 

In order to derive the mean-field F-P equations describing the rotational diffusion, the 
nabla operator — on the left-hand side of Eq. , which is proportional to the linear 
momentum operator and represent the generators of spatial translations, should be replaced 
by the angular momentum operator J representing the generators of rotations: — zV — *■ 
JT" [56], 57 1 . This gives the rotational F-P equation in the following form 

d.f = - [fj, ^ _^ . D„. {/ . .7 . (45) 

where -D^™*'* is the rotational diffusion tensor and the components of the angular momentum 
operator, J, expressed in terms of the Euler angles, u = 7), are given by [ssj l 

( cos (y ^ 

J'l = Jr^ = —i<— COS a cot P da — sin adn -\ \ , (46a) 

I sm 7 J 

J2 = Jy = —i\ — sin a cot (3 da + cos ada ^ \ , (46b) 

I sm 7 J 

J^ = J, = -Ida. (46c) 

When the effective free energy functional is a sum of two term that represent the contri- 
butions coming from the effective internal energy, U[f], and the Boltzmann entropy term 

F[f] = U[f] + {\r^f), (47) 

the variational derivative of the free energy takes the form 

^ = V^ + ln/ + l, y = ^, (48) 

where V is the mean-field potential. 

On substituting the relation ( l48l) into Eq. ( l45l) we obtain the mean-field F-P equation 

dtf = -J ■ D,ot {Jf + f-JV}, (49) 

describing the rotational diffusion governed by the mean-field potential (HHI) . This equation 
can be conveniently cast into the form 

dtf = Dlf + \ [d^UV) + fD^V - VD],f] , (50) 
where the right-hand side is rewritten using the operator 

Dl = -J^D^^J, (51) 



which is quadratic in the components of the angular momentum operator, Jj. 

In the hnear case where the potential, V ^ is independent of the angular distribution 
function, /, the F-P equation fl49l) has been used to study dielectric and Kerr effect relaxation 
of polar liquids based on the rotational diffusion model the rotational motion of molecules 
in the presence of external fields |59|, 0, H H H 

Rotational diffusion of a probe 
molecule dissolved in a liquid crystal phase was investigated in 65,0, 67 . 

When molecules and the orientational distribution function are cylindrically symmetric, 
the model can be described in terms the angle between the electric field and the molecu- 
lar axis lei, whereas angular distributions of a more general form require using both 



azimuthal and polar angles that characterize orientation of the molecules ]6ll . 6^ . In this 
case, for uniaxial (rod-like, calamatic) molecules, the distribution function /(a,/5, 7) = /(w) 
is independent of the Euler angle 7: /(^) = f{4>,6) = /(u), and the angular momentum 
operator can be expressed in terms of the azimuthal and zenithal (polar) angles, and 6, 
as follows 



. ^ dj~^ i L = [r X V] = - [sin Q\ ^e^ 



where 



eg = (r) = cos Q cos x + cos Q sin y ~ sin 6* z, 

= ej^(r) = — sin x + cos y, 
r = 6^(1") = sin Q cos x + sin Q sin y + cos Q z. 



(52) 



(53a) 
(53b) 
(53c) 



When the rotational diffusion tensor is diagonal, D\- — i/j^. 



'i™*'* = dI™*'* and and its elements 
are angular independent, the operator (!^T|) can be written in the simplified form: 

D'j = - [D^r^i^l + D^r^C^l + Df)/:^) . (54) 

In the isotropic case with D^^'^^^ = D^^°^\ we have 



D^j = = D^'°'U [sine]-' deism 9 de) + [sin^]-^^^ 



(55) 



A more complicated biaxial case occurs for asymmetric top molecules |63l | , m acro- 
molecules in liquid solutions [6^ and probes in the biaxial liquid crystal phase [6^. For 
such low symmetry, analytical treatment cannot be simplified and involves the three Euler 
angles, uj = (a,/?, 7). 

Nonlinearity in the lowest order approximation can be introduced through the truncated 
expansion for the internal energy functional U[f] retaining one-particle (linear) and two- 
particle (quadratic) terms 



U[f]= j f/i(a;)/(cj) dcj + i j /(a;i)f/2(a;i, a;2)/('*^2) dcji dLL»2, 



(56) 



da / sin/3d/3 / d7 . . . and f/2(cL'i, (^2) = f^2(<^25<^i) is the sym- 
'0 Jq Jo 
metrized two-particle kernel. The the effective potential 



where / du 



oj[i^> J 



(57) 



is the sum of the external field potential, Ui{uj), and the contribution coming from the 
two-particle intermolecular interactions. 

For rod-like azo-dye molecules, the one-particle part of the effective potential 0571) can be 
written as a sum of the two terms: 

U,{u) = Uiiu) + U,{u), (58) 

where the light-induced contribution 

f//(u) = Uj Eljy ■ Q(u) ■ Euv = Ui luv Qyyiu) (59) 

comes from the interaction of azo-molecules with the activating UV light and the surface- 
induced potential 

Us{u) =UsZ-Q-z = Us Qzz{u) (60) 

takes into account conditions at the bounding surfaces of the azo-dye layer. 
Assuming that the two-particle interaction is of the Maier-Saupe form 

[/2(U1, U2) = U2 Tr[Q(ui) ■ Q(U2)] = U2 Qijiui)Qij{U2) (61) 

we derive the expression for the effective potential of azo-dye molecules 

V{u) = Uj luv <5j/y(u) + Us Qzz{u) + U2Sij Qij{u). (62) 

The equilibrium angular distribution can generally be obtained as a stationary solution 
to the F-P equation (HUI) . It is not difficult to see that the stationary solution given by 

/st(a;) = Z;,' exp[-F(a;)], Z,, = j exp[-F(a;)]da; (63) 

is the Boltzmann distribution determined by the effective potential. Note that the for- 
mula ( l63l) can be obtained from the condition 

where the constant /i can be regarded as a chemical potential that plays the role of the 

Lagrange multiplier defined through the normalization condition j /st(^) da; = 1. 

When the F-P equation is linear, the stationary distribution (163!) describing the equilib- 
rium state is unique. In contrast to the linear case, the effective potential depends on 
the elements of the averaged orientational order parameter tensor ([H]): Viuj) = V{u\S). So, 
the components of the order parameter tensor in the stationary state, Sij = S^j^\ can be 
found from the self-consistency condition 

S^J = J g.,(u)/st(u|S)du. (65) 

In general, there are several solutions of Eq. (1651) representing multiple local extrema (sta- 
tionary points) of the free energy 

F[f St] = i^st(S) = -y S,,S,, - lnZ3t(S). (66) 



Following the line of reasoning presented in Ref. 4J] and using the effective free energy (H7|) 
as the Lyapunov functional, it is not difficult to prove the H theorem for nonlinear F-P 
equations of the form It follows that all transient solutions converge to stationary ones 
in the long time limit. So, each stable stationary distribution is characterized by the basin of 
attraction giving orientational states (angular distributions) that evolve in time approaching 
the stationary distribution. 

Free energy F-P equations (both linear and non-linear) are generally not exactly solvable. 
So, we conclude this section with remarks on numerical methods applicable to nonlinear F-P 
equations. 

The method based on distributed approximating functionals (DAF) which couples the 
path-integral concept to the DAF idea is proposed for numerically solving a general class of 
nonlinear time-dependent Fokker-Planck (F-P) equations in [68^. The approach is applied 
to solve a nonlinear self-consistent dynamic mean-field problem for which both the cumu- 
lant expansion and the scaling theory have been found by Drozdov and Morillo [69| to be 
inadequate to describe a long-lived transient bimodality. 

In Ref. jo^ , a finite-difference method for solving a general class of linear and nonlinear 
F-P equations based on a i^-point Stirling interpolation formula is suggested. A procedure 
to systematically evaluate all the moments of the F-P equation by expanding them in a 
power series in a given function of t is suggested in [t^]. The methods which are extensions 
of this power series expansion formalism to a general Fokker-Planck-Schrodinger process are 
presented in [7l|- They are applied to a well-known problem of the decay of a unstable state 
driven by exponentially correlated Gaussian noise. 



B. Regime of purely in-plane reorientation: 2D model 



In the previous section our model is formulated as the free energy F-P equation fl49p 
describing rotational diffusion of azo-dye molecules governed by the effective mean field 
potential (|62l) . Since general analysis can be rather involved, we first carefully examine our 
model in the limiting two-dimensional case of purely in-plane reorientation. 

When the symmetry is cylindrical and the symmetry axis is directed along the normal to 



the substrates (the z axis), the diffusion coefficients Dx°^'' and Dy°^' are identical: D, 



.(rot) 



D 



(rot) 



In this case, the expression for the operator (13^ is given by 

dl\+D^r''^dl c = cos^. 



D 



(rot) 



1-C2 



where 



The simplified two-dimensional model 

d^f =d^[d^f + fd^V 

+ fdlv - vdlf 



c')d,] = {l-c')d',-2cd, 



dlf + l[dl{fv) 



(67) 



(68) 



(69) 



immediately follows from Eq. fl^Uj) provided that D^" = and /(u, t) = f{(l),t). 



More interestingly, the F-P equation (lUUl) can be derived by assuming that the out-of- 
plane component of the unit vector u describing orientation of the azo-dye molecules is 
suppressed and, as is shown in Fig. [H u = (sin0, cos0, 0). It implies that the molecules 
are constrained to be parallel to the substrate plane (the x-y plane) and the orientational 
distribution function takes the factorized form 

/(u,t) = /(0,t)5(c), c = cos^, (70) 

where 5{c) is the 5-function. 

Derivation procedure involves two steps: (a) substituting the relations (!70l) and (1671) 
into Eq. ( |50i) ; and (b) integrating the result over the out-of-plane variable c. [Recall that 



du = 
the relation 



27r 



1 y dc] . The kinetic equation (!69l) can now be readily recovered by using 



dc = 0. 



(71) 



So, in this section we consider the F-P equation (1691) representing the simplest case 
when the out-of-plane reorientation has been completely suppressed. Then, for u = 
(sin(^, cos(/), 0), the diagonal elements of the order parameter tensor ([T]) averaged over the 
azimuthal angle are given by 

1 





{Qxx) 


Sy = 


{Qyy) 







;(3sinV- 1) 

■(3cosV- 1) 
1 

2' 



-3(cos2</)) + 1 
3 (cos 20) + 1 



(72a) 
(72b) 
(72c) 



2tt 



where (...) = 

We can now substitute the order parameters (1721) into the effective potential V given in 
Eq. The result for the angular dependent part of the potential is 

V = {vi + f 2 (cos 20) ) cos 2(f) = V cos 20, (73) 

where Vi = Sujluv/^: and V2 = 9^2/8 are the light induced and intermolecular interaction 
parameters, respectively. 

Our next step is to obtain the system of equation for the averaged harmonics, c„(r) = 
(cos n0)(r), that are proportional to the Fourier coefficients of the distribution function, 
/(0, r). To this end, we integrate the F-P equation (!69l) multiplied by cos ruj) over the 
azimuthal angle and apply the relations 

f.27r 

COS?T,09|/d0 = ((9^COSn0) = —n^ (cos 77,0) = — n^Cn, 

r.27r 

COS 77,0 5|(/y) d0 = (Vd'^cosncp) = — in? (y cos ncj)) = —n^v (c„+2 + Cn-2)/2, 

r.27r 



COS n(j)fd^Vd(f)= (cos n0 d^V) = -2v (c„+2 + c„_2) , 



2tt 



COS 710 V d^f I 



{dl{V cosn, 



--[{n + 2fcn+2 + {n-2fcn-2]- 



(74a) 
(74b) 
(74c) 
(74d) 



The resulting system reads 

drCn = -n'^Cn + nv{Cn+2- Cln^2\), Cq = 1, (75) 

where v = Vi + V2C2- 

When /(0 + tt) = the odd numbered harmonics vanish, C2k+i = 0. For the even 

numbered harmonics, pk = C2k, the system ( |75l) can be conveniently recast into the form 

drPk = -4:k'^Pk + 2kv{pk+i - Pk-i), k = 1,2,... (76) 
Pk = C2k, Po = Co = 1, V = vi+ V2P1, (77) 

where pi = (cos 2(f)) is the order parameter harmonics that enter the expressions for the 
orientational order parameters fl7^ . 
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FIG. 2: (a) Intersection points represent solutions of self-consistency equation (18^ . 
(b) Stationary state free energy (166|1 as a function of the parameter v. 



1. Bifurcations of stationary states 
From the general formula (1631) we obtain the expression for the stationary distributions 

l'2n 

fst = Z-^exp[-V]= Z;,'exp[-vcos2(l)], Z,,= exp[-r]d0 (78) 

Jo 

representing the photosteady states in the two dimensional case with the potential (1731) . 
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FIG. 3: (a) Order parameter pi = (cos 20) as a function of the parameter vi at V2 = —2.5; 
(b) bifurcation curves in the f2-f i plane are typical of the cusp catastrophe with the cusp 

singularity located at (—2,0). 



Equation ( 1781) can now be combined with the relation 

oo 

exp[-t;cos20] = /o(t;) + 2 ^(-l)'=4(t;) cos2A;0, 



(79) 



k=l 



where Ik is the modified Bessel function of integer order j72|, to derive the formulas 

Zst = 27r/o(t;), (80) 
= {-irh{v)/Io{v). (81) 

giving the stationary state statistical integral, Zst and the averaged harmonics, p^^*'*, ex- 
pressed in terms of the parameter v. Using the recurrence relation [7^] 



v[Ik^i{v) - Ik+i{v)] = 2klk{v) 



(82) 



it is not difficult to verify that the formula (!HT|) gives the stationary solution to the sys- 
tem (TfGl) which, in the steady state regime with drPk = 0, is represented by the finite 
difference equation 



(83) 




FIG. 4: Bifurcation diagram as the cusp surface in the three dimensional (fi,f2,Pi) space. 

From Eq. ( 1731) . the parameter v = Vi +f2pf*'' depends on the order parameter harmonics 
and Eq. ( ISTl) with k = 1 provides the self-consistency condition 

) = {v- vi)/v2 = -h{v)/Io{v). (84) 

This condition can also be obtained as the stationary point equation for the stationary state 
free energy (IMj) . In our case, we have 

f/[/st]=-#^-|bn'^ (85) 

{lnU) = -lnZ,,-vp^f\ (86) 

Fst{v) = v^'[-^ + v,v] -\nIo{v), (87) 

where the additive constant is chosen so as to have the free energy vanishing at w = 0. 

In Fig. [2^a) it is illustrated that, in the v-pi plane, solutions of the self-consistency 
equation flH^ can be found as intersection points of the curve —Iiiy) / I^iy) and the straight 



line, {v — Vi)/v2. It is seen that the number of the intersection points varies between one 
and three depending on the values of the parameters Vi and V2. 

As is shown in Fig. Mjo), for the case of three stationary states, the free energy curves 
are of the double-well potential form with two local minima separated by the energy barrier. 
From the lowest order term of the series expansion 

-/i(x)//o(x)^ -- + --- (88) 

it is not difficult to see that this case may occur only if the parameter V2 is less than —2. 

Referring to Fig. [21(b), at f i = and V2 < —2, the free energy (lH7j) is an even function 
of V, Fst{v) = Fst{—v), with two symmetrically arranged minima representing two stable 
stationary states. When the parameter of intermolecular interaction, V2, increases passing 

(c) 

through its critical value, V2 = —2, the minima come close together and coalesce at the 
critical point. So, at V2 > —2, there is only one minimum corresponding to the unique 
equilibrium state. 

When the activating light is switched on, the parameter vi is distinct from zero. It gives 
rise to asymmetry effects illustrated in Fig. [21(b). It can be seen that, at Vi 7^ 0, one of two 
minima becomes metastable. The local maximum representing the unstable stationary state 
and the metastable minima merge and disappear provided the magnitude of the parameter 
vi is sufficiently large. 

This effect is evident from the curve depicted in Fig.[3]^a) where the stationary state order 
parameter harmonics is plotted in the vi-pi plane by using the following parametrization 

„(st)_U=Pl(0 = -/l(0//o(0, 

~U=vm=^-V2M0. ^ ^ 

So, the free energy has two local minima only if the inequalities 

V2 < = -2, f < fi < (90) 

are satisfied. The critical values of the parameter vi depend on the reduced strength of 
intermolecular interaction f 2 and can be parameterized as follows 



v2 = v2{o = -[{h{o/my,r' 



(91) 



i+Pi(0/e-p?(0' 



Geometrically, in the V2-V1 plane, equation (]9Tl) defines the bifurcation curves shown in 
Fig. [31(b). These curves form a bifurcation set which is the projection of the cusp surface 

Sb=Iv2 = C (92) 

[pi = -hio/m 

representing the bifurcation diagram in the three dimensional {vi,V2,pi) space (see Fig.Hj). 
Note that the cusp bifurcation occurs as a canonical model of a codimension 2 singularity [73| 



and the surface shown in Fig. [H is typical of the cusp catastrophe [T^, [75 



We conclude this section with the remark on how diffusion models may account for 
the effect of long-term stability by using approximation of the "frozen" potential proposed 
in [43^1. Mathematically, it implies that, after switching off the exciting light at time t = t^s 
with the order parameter harmonics Pos = Pi(tofr), the relaxation process is governed by the 
kinetic equations for the harmonics (|76l) where the parameter v is changed to the "frozen" 
interaction parameter Vos = V2PoS- From Eq. (IHTl) the stationary value of the order parameter 
harmonics 

pis^ = -hM/IoM (93) 

is determined by the "frozen" strength of intermolecular interaction, Vos, and thus depends 
on the value of the order parameter harmonics at the time of switching, t = tos- So, in the 
two state and in 2D diffusion models the memory effect underlying the long-term stability 
of POA is described by the relations (j38l) and (l93ll . respectively. 



IV. RESULTS 

In Sec. Ill CI and Sec. IIIIB| we employed the approaches based on the master and Fokker- 
Planck equations to introduce two different models: the two state model and the two dimen- 
sional diffusion model, respectively. In both cases, the photoinduced anisotropy is charac- 
terized by the orientational order parameters whose temporal evolution is governed by the 
kinetic equations of the model. 

In Sec. Ill Al we discussed how the order parameters can be related to absorption char- 
acteristics such as extinction (absorption) coefficients and optical densities. Specifically, 
Eq. (fT3l) shows that the absorption order parameters, Sl"'\ defined in Eq. ([8]) as a function 
of the principal values of the extinction coefficients, k,, are proportional to the orientational 
order parameters (0). So, a comparison between the theory and experiment can be made 
from measured values of the absorption coefficients. 

In thin anisotropic films, the absorption coefficients can be determined experimentally 



using the methods of ellipsometry |76l. l77l|. These are generally based on the analysis of the 
polarization state of light reflected from or transmitted through a sample. 

One of the simplest experimental procedures is to measure the light transmittance of 
a film when the testing beam is normally incident and linearly polarized. Performing the 
measurements for beams polarized perpendicular and parallel to the polarization vector of 
the UV light the two in-plane optical densities, Dy"'^ = D^^^^ and D^^ = D^f^ , can be obtained 
function of the irradiation dose. 

The normal component, Di"'\ then can be estimated by assuming that the total sum of 
principal optical densities 

Dl:l = D^^ + D^^ + Df^ (94) 

does not depend on the irradiation dose and the photosaturated state is uniaxially anisotropic 
with D^^ = D^f^ . More details about this approach can be found, e.g., in Refs where 
it was applied to azopolymer films. 

In the Appendix we show that the absorption extinction coefficients can be extracted 
from the dependence of absorbance on the incidence angle measured using probe beams 
which are linearly polarized parallel (p-polarization) and perpendicular (s-polarization) to 
the plane of incidence. In order to fit the experimentally measured curves, this method 
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FIG. 5: Extinction coefficients as a function of irradiation time. 



relies on the analytical expressions for the transmission coefficients of biaxially anisotropic 
absorbing layers deduced in the Appendix (see Eq. flA28p ). 

The results for the extinction coefficients, Kx, K,y and summarized in Figure 

where the coefficients are plotted against the irradiation time. The corresponding absorption 
order parameters, Si^\ Sy""^ and Si"'\ evaluated from the experimental data by using the 
formula (IHD are presented in Fig. [6] and Fig. [71 

It can be seen that the initial and photosaturated states are both uniaxially anisotropic 
with Si""^ = Sy'^'' and Si""^ = Sy"'\ respectively. So, similar to the case of azopolymers, the 
transient photoinduced orientational structures are inevitably biaxial. 

It is also clear that, before reaching the regime of photosaturation, the in-plane order 
parameters, Si""^ and Sy"'\ undergo pronounced changes. By contrast, the normal component 
of the order parameter, Si"'\ slowly increases with irradiation time. We can therefore employ 
the two dimensional diffusion model described in Sec. IIII Bl as a zero order approximation 
where variations of the normal order parameter component, Si"'\ are neglected. 

The theoretical curves shown in Figure as solid lines are computed by solving the 
system flTHj) . The ffiting procedure is as follows. 

Assuming that the order parameter S"!"^ is constant and Si"'"' ~ —0.334, we obtain the 
coefficient that enter the relation (I13p . ^ 0.67, and the absorption anisotropy parameter 
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FIG. 6: Absorption order parameters as a function of irradiation time. The theoretical 
curves are computed by solving the system of equations ( !76l) for the two dimensional 

diffusion model. 



cTq ~ 6.1. Then, from the experimental data, we can estimate the order parameter harmonics 
in the photosteady state, pf*'* ~ —0.887 ( the corresponding in-plane order parameter is 
Sxx^ ~ 0.9 ). Substituting this value into the self-consistency condition flH^ gives the 
equation linking the two dimensionless interaction parameters: Vi and V2- At fi = 1.0, 
this equation can be solved to yield the value of the intermolecular interaction parameter, 
V2 ^ -4.22. 

According to the experimental data presented in Figs. [5]- [71 the irradiation time it takes to 
reach the regime of photosaturation is about 64 min. From the other hand, for the computed 
dependence of the order parameter harmonics, pi, on r = D^r'h, this regime takes place at 
r > 1.1. So, the rotational diffusion constant Z)i™*'* can be estimated at about 0.017 min~^ 

2.8 X 10"^ s-i). 

Referring to Fig. [HI agreement between the theoretical curves and experiment indicates 
that the two dimensional diffusion model can be regarded as a good approximation to start 
from. So, the regime of kinetics of the photoinduced structures in the azo-dye film appears 
to be close to the limiting case of the in-plane reorientation. 

Now we consider the two state model formulated in Sec. Ill CI Similar to the case of 
the 2D diffusion model, our first step is to determine the coefficient ra and the anisotropy 
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FIG. 7: Absorption order parameters as a function of irradiation time. The theoretical 
curves are computed by numerically solving the kinetic equations for the two state 

model (EI]). 



parameter aa- The coefficient r^, ^ 0.89, can be calculated as the solution of the equation 
obtained by substituting the photosaturated value of the the absorption order parameter, 
SifJ = TaSst ~ 0.59, into Eq. (H^ . We also find that the absorption anisotropy parameter 
aa is about 25.1, cTq ~ 25.1, and the photosteady state is characterized by by the order 
parameter Sst = Sxx^ ~ 0.66. 

The numerical results shown in Fig. [7] are computed in the regime of photoreorientation 
where the decay rate of the excited state 7e ~ 2.5 min~^ is much larger than the excitation 
rate qgluv/le ~ 10"^, and the thermal relaxation rates: 7G/7e ~ 5 x 10~^ and ^mlle ~ 
2 X 10~^. Numerical calculations in the presence of irradiation were followed by computing 
the stationary values of S and AS* to which the order parameters decay after switching off 
the irradiation at time tofr (for more details see Sec. Ill C ip . 

The results presented in Fig. [7] suggest that the two-state model can be used to relax the 
assumption on purely in-plane photoreorientation and to go beyond the limitations of the 
2D diffusion model. 



V. DISCUSSION AND CONCLUSIONS 



In order to study the kinetics of the photoinduced ordering in azo-dye films we employed 
two different models. 

The two state model is formulated by using the phenomenological approach developed in 



Refs. [35|, |40[ . In this approach, the film is represented by an ensemble of two-level molecular 
systems. So, it starts from the master equation ( ITTI) for one-particle angular distribution 
functions of the ground-state and excited molecules. The kinetics is then determined by the 
angular redistribution probabilities that enter the photoexcitation and decay rates. They 
also define coupling between the azo-dye molecules and the anisotropic field represented by 
the distribution function of the matrix fm- This anisotropic field reflects the presence of 
long-living angular correlations and stabilizes the photoinduced anisotropy. 

The resulting kinetic equations ( |3T]) for the order parameter components are deduced by 
using the parabolic approximation suggested in Ref. [35| to express the order parameter cor- 
relation functions in terms of the order parameter tensor. Following the procedure described 
in Refs. [Isl, 3^, these order parameter correlation functions are additionally modified in 



order to take into account constraints suppressing out-of-plane reorientation. Another im- 
portant assumption taken in our two state model is that the excited molecules are isotropic 
and do not affect the ordering kinetics directly. 

Similarly, there is an alternative approach which is formulated in Sec. IIIII without ex- 
plicit reference to excited electronic levels. According to this approach, the photoinduced 
anisotropy arises from the rotational Brownian motion of azo-dye molecules in the effective 
light modified potential. Mathematically, this suggests using the mean-field Fokker-Planck 
equation ( H5i) with the effective free energy functional (H7I) as the equations governing the 
kinetics of photoinduced ordering. Thus diffusion models can be defined by specifying the 
rotational diffusion tensor and the effective potential (jUj) that enter Eq. (1491) . 

The two dimensional model studied in Sec. IIII Bl presents the simplest case to start from. 
It is based on the approximation of purely in-plane photoreorientation which assumes the 
normal order parameter component Sz kept constant. 

In order to test applicability of this approximation, we compared the predictions of this 
simple model with the available experimental data. Fro Fig. [6] it is clear that, in azo-dye 
films, the kinetics of photoinduced structures take place in the regime close to the limiting 
case of purely in-plane photoreorientation. 

Referring to Fig. [6] and Fig. [71 the comparison between the numerical results and the 
experimental data shows that the two-state and the 2D diffusion models both correctly 
capture the basic features of the photo-ordering kinetics in the azo-dye layers. It comes as 
no surprise that the results computed from the two-level model give better agreement with 
experiment than the ones for the 2D diffusion model. The primary reason for this is that 
the two-state model takes into account effects due to variations of Sz- 

These effects can also be taken into consideration in the rotational diffusion approach by 
expanding the orientational distribution function into a series over the spherical harmonics. 



Yimiu) or, more generally, the Wigner D functions [58|, 178[, D^^,(a;). The mean-field 
Fokker-Planck equation (150|) then can be transformed into the system of nonlinear ordinary 
differential equations for the averaged harmonics, {Yim)(t) or (-D^^,)(t). Equations (ITHjl 
represent the special case of such system derived for the two dimensional model. 

For an infinite number of equations, numerical analysis involves truncating the system 
so that only a finite number of harmonics are taken into account. The number of harmon- 



ics is typically determined by the required accuracy of calculations. Difficulties emerge if 
this number turns out to be very large. For instance, this is the case for highly ordered 
photosteady states. 

So, we have demonstrated that the phenomenological approach of Ref. and gen- 
erahzed diffusion models can be used as useful tools for studying photoinduced ordering 
processes in azo-dye ffims. It should be noted, however, that theoretical approaches of this 
sort, by definition, do not involve explicit considerations of microscopic details of azo-dye 
film physics. A more comprehensive study is required to relate the effective parameters of 
the models and physical parameters characterizing interactions between molecular units of 
films. 
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APPENDIX A: TRANSMISSION COEFFICIENTS OF BIAXIALLY 
ANISOTROPIC ABSORBING LAYERS 

In this section we derive the exact solution to the transmission boundary value problem 
by applying the theoretical approach developed in Refs. 0, [si to the case of biaxially (and 
uniformly) anisotropic absorbing layers. 

As is shown in Fig. [8], we consider an absorbing uniformly anisotropic film of thickness 
d with the z axis giving the optic axis normal to the bounding surfaces: 2; = and z = d. 
The other two in-plane optic axes are assumed to be directed along the unit vectors x and 
y. In this case the dielectric tensor of the ffim is diagonal and is defined in Eq. (IH]). From 
Eq. (171), the principal values of the tensor, ea, can be expressed in terms of the refractive 
indices, na \ and the extinction coefficients, Ka- 

The medium surrounding the layer is assumed to be optically isotropic and characterized 
by the dielectric constant Cm and the magnetic permittivity /im- Referring to Fig. [HI there 
are two plane waves in the half space z < bounded by the entrance face of the layer: the 
incoming incident wave {Einc, Hinc} and the outgoing reflected wave {Erefi, Hrefi}. In the half 
space z > d after the exit face, the only wave is the transmitted plane wave {Etrans, Utrans} 
which propagates along the direction of incidence and is excited by the incident light. 

So, the electric field outside the layer is a superposition of the plane waves 

EUo = Ei,e(kinc)e'^'^^'-'-) + E,efl(k,efl)e^(^-«-'-), (Ala) 
E|,>, = Etr(ktr)e^^''-"-), (Alb) 

where the wave vectors kinc, kj-efi and k^i- that are constrained to lie in the plane of incidence 
due to the boundary conditions requiring the tangential components of the electric and 
magnetic fields to be continuous at the boundary surfaces. These conditions are given by 

z X [E|^=o+o - E|^=o-o] = z X [B\,=d+o - ^z=d-o\ = 0, (A2a) 
z X [H 1^=0+0 - H|^=o-o] = z X [ii\z=d+o - ii\z=d-o\ = 0. (A2b) 



Another consequence of the boundary conditions flA2p is that the tangential components 
of the wave vectors are the same. Assuming that the incidence plane is the x-z plane we 




FIG. 8: Geometry of anisotropic layer in the plane of incidence. 



have 

ka = A;vacqa = k^J^a = Ki + fc^"^ z, a E {inc, refl, tr}, (A3) 

where km/Kac = = y/fim^m is the refractive index of the ambient medium and fcvac = uj/c 
is the free-space wave number. The wave vector components can now be expressed in terms 
of the incidence angle 6'inc as follows 

kx = sin ^inc = ^vacQ'x; (^4) 
kz = ^ = k[, ^ = fcm cos ^inc = ^vac Q'm; (-^5) 

Qx = sin 6'inc, = - Qx- (A6) 

The plane wave traveling in the isotropic ambient medium along the wave vector (1A3P is 
transverse, so that the polarization vector is given by 

E„(k,) = Ej"^ei(k„) + E'f^e^ika), (A7) 
ei(k,) = A;-i(A;(")x-fc,z), e2(k„) = Y, (A8) 





FIG. 9: Absorbance as a function of the incidence angle. Circles and squares represent the 
data measured for testing beam linearly polarized parallel (p-polarized wave) and 
perpendicular (s-polarized wave) to the plane of incidence, respectively. The theoretical 

curves are shown as solid lines. 



where £"11 (= Ep ) and E " (= Es ) are the in-plane and out-of-plane components of the 
electric field, respectively. The vector characterizing the magnetic field is 

fi^ H,(k,) = X E„(k„) = [^{"^y - E^f'^e.iK)] , (A9) 

where = k^^J^a = n^^a- Note that, for plane waves, the dimensionless vector 

q = K^k (AlO) 

is parallel to k and its length gives the refractive index. For convenience, we shall use this 
vector in place of the wave vector. 

The electromagnetic field of incident, transmitted and reflected waves propagating in the 
ambient medium is of the general form 

{E, H} = {E{z), H{z)} e*('=-^-'^*). (All) 



On substituting the relations (lAlip into the Maxwell equations we can obtain the equations 
for the tangential components of the electromagnetic field inside the anisotropic layer. The 
result can be written in the following 4x4 matrix form 

-,3.F^M.F.(^^^'A(^'-], ....... (A12) 



where Ep(.) = j;)) and Hp(.-) = (_^(-;^)) 



For the dielectric tensor ([6]) with the plane of incidence parallel to the x-z plane, from 
the general expressions derived in Refs. [79|, l80| , the 2x2 matrices Mjj characterizing the 
block structure of the matrix M are given by 

Mn = ej' ,"2) , M„ = 0, (A13) 

M., = ("o" , M.. = 0. (A14) 



According to the computational procedure developed in Refs. [79|, ISOJ , the transmission 
and reflection matrices defined through the linear input-output relations 




^(inc)' 



(A15) 



s 



R I ^!inc) I (A16) 



can be expressed in terms of the linking matrix 



w = v-.u-'w.v„ = (^;;^-l (A17) 



as follows 



T = W^,\ (A18) 
R = W21 ■ W^-^^ = W21 ■ T. (A19) 

The expression for the linking matrix ( 1A17I) involves the inverse of the evolution operator 

V~\h) = V{-h) = exp{-zM h}, h = K^d (A20) 

and the eigenvector matrix for the ambient medium 



V. = T^^" (A21) 



Em — 

which is characterized by the two diagonal 2x2 matrices 

Em = diag(gm/nm, 1), /^m Hm = diag(nm, ^m), (A22) 



where 0-3 = diag(l, —1). 

In our case, the resulting expression for the evolution operator is 



U(/^)=exp{^MM = V.(^^+ G=(^j°^{jj, (A23) 

where I2 = diag(l, 1) and 

U± = exp{±2Q h}, Q = (^^^ °J , (A24) 



H^i(-°). ,A2T) 

We can now substitute the operator (1A23I) into the linking matrix (1A17I) and obtain the 
transmission and reflection matrices using the relations (1A18I1 and (]A19[) . The result is given 

by 



tpiqx) \ _ I 



ts{q^)J l2-U^P 



-Qm ± Qp 



V-V;\ V±=hur2^"'" /i^n,"^ . (A30) 



/i ± /im^gmy 



From Eq. flA28p and Eq. flA29p . non-diagonal elements of both transmission and reflection 
matrices vanish. Algebraically, this is a consequence of the diagonal form of the block 
matrices that enter the operator of evolution flA23p . 

So, absorption of plane waves linearly polarized parallel and perpendicular to the plane 
of incidence can be characterized by the effective optical densities. Dp and Dg, expressed in 
terms of the corresponding transmission coefficients: 

Dp,sie,nc) = Dp^siQx) = -21n|tp,,(g,.)|- (A31) 

The optical densities ( ]A3ip are proportional to the absorbances measured experimentally, 
Dp^^''^ and Di'^^'^\ and determine the theoretical dependence of the absorbance on the inci- 
dence angle, 6'inc (or, equivalently, on the incidence angle parameter = n^smOinc)- 

In Fig. [9], the experimental data on angular dependence of absorbance measured in the 
azo-dye SDl film of the thickness 15 nm at different irradiation doses are fitted by the 
theoretical curves computed from the formula flA3ip . 



For our purposes, full description of a rather standard experimental procedure is not 
important (more details can be found in [12]). So, without going into details we note that 
the film was illuminated with linearly polarized UV light at varying exposure time by using 



LED exposure light source. The wavelength and the intensity of the actinic light were 
365 nm and 3.0 mW/cm^, respectively. 

In Fig. El the extinction coefficients of the azo-dye layer found as the ffiting parameters 
are plotted as a function of the irradiation time. 
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